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' Abstract The Magnificent Seven have all been discov- 
' ered by their exceptional soft X-ray spectra and high ra- 
, tios of X-ray to optical flux. They all are considered to be 
nearby sources. Searching for similar objects with larger 
distances, one expects larger interstellar absorption re- 
sulting in harder X-ray counterparts. Current interstel- 
lar absorption treatment depends on chosen abundances 
and scattering cross-sections of the elements as well as on 
the 3D distribution of the interstellar medium. After a 
discussion of these factors we use the c omprehensive 3D 
measu rements of the Local Bubble by ILallement et al.l 
Hods') to construct two simple models of the 3D distri- 
bution of the hydrogen column density. We test these 
models by using a set of soft X-ray sources with known 
distances. Finally, we discuss possible applications for 
distance estimations and population synthesis studies. 
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1 Introduction 

The Magnificent Seven (M7), as the ROSAT-discovered 
X-ray thermal isolated neutron stars are sometimes dubbed 
are exceptional because of their soft blackbody-like spec- 
tra without non-thermal com ponents; for r eviews see e.g. 
Haberl 2006 (this volume) or Haberll l|2004|) . One of them 
has the second-best blackbody spectrum after the cosmic 
background radiation and this lead to o ne of the best 
neutron star radius estimates known f^Triimper et alJ2004l) 
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Thought to be nearby, the M7 represent nearly half of 
the local young neutron star population. Despite sev- 
eral searches for n ew candidates of this radio-quiet neu- 
tron star class (e.g.'Agrieros et al.l2006tfChieregato et al.l 
12005: Rutlcdge et al. 2003), no new M7-like object has 
been confirmed up to now. One aggravating circumstance 
these searches face is the hardening of the X-ray spec- 
trum due to the interstellar absorption. Caused mainly 
by photo-electric absorption of photons by heavy ele- 
ments, the interstellar absorption is usually described 
by the equivalent hydrogen column density assuming cer- 
tain chemical abundances. Highly important is the dumpi- 
ness of the interstellar medium (ISM) causing (un) favourable 
lines of sights to search for new M7-like objects. The 
knowledge about the local ISM structure is also impor- 
tant in connection with the long-debated issue of isolated 
neutron stars accreting from the ISM. 

1.1 Abundances and cross-sections 

Most absorption models today presume element abun- 
da nces independent of the line of sight. It was outlined 
bv lWilms et all lj2000(l that the total gas plus dust ISM 
abundances are lower than the local - Solar - abundances 
and the general ISM abundance uncertainties are still in 
'the order of 0.1 dex because the measurements are very 
difficult. One can choose between several different abun- 
dance tables implemented e.g. in XSPEC for analysing 
X-ray spectra. We note here only t wo examples com 



moniy used - the abundance tables by^bide^^^Grevess^ 
1119891) ('angr' in XSPEC) and bv lWilms et all i200nD 

('w ilm' in XSPE C) considering newer measurements by 

e.g.lSnow fc Witt (1996) ; CardeUi et aL (19M) ; Mever et al 
lll997lll99SD . 

The individual abundances and photoionization cross- 
sections of the elements and their ions as well as dust 
grain properties (e.g. sizes, shapes) influence the total 
photo-electric absorption of X-rays. Again, we mention 
here only two examplar y works used very often for anal- 
ysis of X-ray spectra. iMorrison fc McCammonI (|l_983') 
did a polynomial fit of the effective absorption cross- 
section per hydrogen atom based on measured atomic ab- 
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Fig. 1 The effect of different abundances and cross- 
sections on the absorption towards soft X-ray sources 

For the same neutron star (blackbody of 90 eV, 11.67 km ra- 
dius, and 1.48 Solar masses) the ratio of the X-ray fluxes (as 
would be expected to be observed by ROSAT) resulting from 
the two XSPEC routines Hbabs' and 'wabs' are calculated 
for different distances. The galactic coordinate x varies from 
50 to 2500 pc (direction: from the Sun towards the galactic 
center), y is always 0, z is 25 pc north . For the ISM distri- 
bution, we applied the same as used bv lPopov etldj 11200011 . 
which is described in more detail below. The effect of less ab- 
sorption, resulting from the use of the newer routine 'tbabs', 
is stronger at larger column densities, hence distances. 




sorption cross-secti ons bvlHenke et all lll982ll and abun- 
dances mostly from [^ide^^^^^^^^_^^^) . With the 
exception of oxygen jMorrison fc McCammonl l|l98t^ as- 
sume the elements to be either entirely in the gas phase 
or completely depleted into du st grains. The recen t im- 
proved absorption treatment bv lWilms et alJ l)2000() takes 
into account newer element abundances as noted above 
and more r ecent photoionizati o n cross-section calc ula- 



Fig. 2 The limits of using the lFerlet et alj II1985I) for- 
mula 

T he sodium column densi ty as derived from the density data 
by ILallement et al] (1200311 is shown in galactic cartesian co- 
ordinates at z=0; x is towards the galactic center which is at 
the right, y towards 1=90° is pointing up. The image spans 
500 pc times 500 pc. The area within the cont our has low 
column densities (log N(Nal ) < 11) for which Weltv et al.l 
(1994) doubted whether the IFerlet etaT. (1985) formula is 
applicable. This region of small sodium column densities is 
largest at z=0 and becomes smaller for lower or larger z. 



tions by e.g. 



IWilms et al 



Verner et al1(|l993j) : l\^ner fc YakovlevI (( 199 



ular cross-section f or H2. Besides this updated database 
IWilms et al.] (1200(1) considered furthermore a simple spher- 
ical composite dust grain model in an MRN distribution 
(|Mathis et alJll977|) . 

Both models are implemented in XSPEC and can be 
used via the ro utines ^wabs^ (iMorrison fc McCammonl 
Il983|) or ' tbabs' (|Wilms et alJl200(]|) . The various absorp- 
tion descriptions can lead to different results (see Fig.^l 
or expectations of e.g. observable objects. The differences 
caused by the diverse treatments show our incomplete 
knowledge about the interstellar absorption and should 
be kept in mind especially when interpreting soft X-ray 
spectra. 

1.2 The inhomogeneously distributed ISM 

In X-ray astronomy with low spectroscopic resolution 
the chosen abundances and cross-sections are used to 
derive the hydrogen column density responsible for the 
instellar absorption. Cold and warm H-atoms (Hi ) as 
well as hydrogen molecules and ionized hydrogen (Hll 
) contribute to the overall amount of hydrogen along a 
line of sight. These different components of the hydrogen 
column density are estimated by different measurement 



methods, e.g. various extinction measurements as sum- 
marized by .Knu(ie (2,0 02.) or radio 2 1-cm observations of 
hydro gen as the extensive stu dy bv iDickev fc Lockmai] 
||199(]|) and more recently by iKalberla et all (|2005|) . A 
wealth of theoretical models has been applied to these 
(2000) additionally include an improved molec-°b™tions. However, most of th em are only 2D with 

tew more recent ex ceptio ns (see e.g.lDrimniel et al.l2003l 
lAmores fc Lepind l2005l or iMarshall et all l2006l for~3D 
models). The overall result is a highly inhomogeneously 
distributed interstellar medium within the Milky Way, 
present in bubbles with loops or dense rims, and shaped 
by molecular clouds, supernova explosions and stellar ac- 
tivity (for an illustrated overview seeiHcnbcst & CoupgJ 
Il994|) . The ISM distribution is best known in the close 
Solar neighbourhood, especially in the Local Bubble (e.g. 
^rcitschwcrdt 1998). The most detailed view is provided 



by the study of La llement et all ^200^ . who measured 
Na I absorption towards 1005 sightlines with Hipparcos 
distances up to 350 pc from the Sun; for more details 
on the measurements and the applied density inversion 
method see Lallement et aL L2j)Q3) , the prec ursor paper 
ISfeir et al.1 |l999|l . and I Vei^elv et all (|200l[l . Nai is re- 
garded as good tracer of the total amount of the neu- 
tral interstellar gas and the sodium column density is 
thought to be convertable in to a hydrogen column den- 
sity (e.g. iFerlet et alJ () 1985ft . but see also discussion in 
sec. O. For more distant regions the determination of 
a 3D distribution of the ISM is hampered mainly by 
the unknown or imprecise distances of the measured ob- 
jects and the lack of a well-sampled grid of measurement 



The Magnificent Seven in the dusty prairie 



3 




+ 90 deg Galactic longitude + 270 deg 



Fig. 3 The N(H) predicted by the analytical model 
at 1 kpc 

Positive galactic latitudes are up and / = 0° is in the center. 
The N(H)is in units of lO^^cm"^. For description of the model 
and discussion see text. 
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Fig. 4 The N(H) predicted by the extinction model 
based on Hakkila et al. (_1997) at 1 kpc 

Positive galactic latitudes are up and Z = 0° is in the center. 
The N(H) scale is in unities of 10^^cm~^. For description of 
the model and discussion see text. 



points. This is true for all extinction-related methods 
to our knowledge and results in much more coarse sam- 
pling compared to the Solar neighbourhood. The results 
obtained by the various methods can differ for the same 
regions and a combination can be difficult. These diffi- 
culties result partly from the different observed astro- 
physical phenomena like measuring the dust emission or 
stellar magnitudes. They also result from different sets of 
stellar absolute magnitude calibrations and the underly- 
ing assumptions which can deviate from each other as in- 
terstellar extinction determination is an iterative process 
coupled with the identification o f the statistical proper- 
ties of stars (iHakkila et al.lll997|) . 



2 Two simple 3D distribution model cubes of 
the absorbing ISM 

Considering t he absorption in the So lar neighbourhood, 
the result bv iLallement et all l)2003(l is a good starting 
point representing the currently best database of the lo- 
cal ISM distribution. R. Lallcment kindly provided us 
with the Nal density cube derived using the inversion 
method developed bv IVergelv et al.l l)200l|) . Due to the 
smoothing length of 2 5 pc applied in this method by 
ILallement et al.l l)2003() . structures smaller than 25 pc 
cannot be resolved. We note further that even if measure- 
ments go up to Hipparcos distances of 350 pc from the 
Sun, the sampling becomes coarser at larger distances. 
Thus the sodium density cube has a span of only 250 pc 
and starting from 200 pc one has to be careful dealing 
eventually with the a priori density information applied 
in the inversion method (R. Lallement, personal com- 
munication). From the Nai density cube we calculate the 
column density for the same grid (sampling 3.9 pc). This 
is then c onverted to H i colu mn density applying the for- 
mula bv lFerlet etall l)l985D . 

It has to be mentioned that there are on-going discus- 
sions about how well the sodium D-line absorption ac- 



tually traces H i . The correlation found between the H I 
and Na column densities d erived from Na- D absorption 
lines b v lFerlet et aP l)l985 ') was doubted bv lWeltv et alJ 
(|l994l) . especially for low column densities (log N{NaI) < 
11). The region of the Lallement cube that should not 
be used because it has such low column densities is indi- 
cated in Fig. 121 

However, the lowest known X-ray-measured hydro- 
gen column density for one of the M7 converts to a 
sodium column density of log N(N ai) > 11.5 (applying 
the formula by iFerlet et alJ Il985() . Therefore this low- 
column density uncertainty is not important consider- 
ing the M7 and negligible for the intended population 
synthesis with neutron stars (NSs) at usually larger dis- 
tances^^ 

IVergelv et al.l l)200ll) could not always find a correlation 
of the Nai density with the H i density. They concluded 
that the correlation is rather weak due to non-constant 
population ratios for interstellar med i tim sp ecies, thus 
different abundances. ILallement et alJ l)2003|) noted the 
size of the Local Bubble revealed by H i is smaller than 
when derived by Nai. This was explained by the first ion- 
ization stage of Nai being below tha t of H i , resultin g 
in a long er neutral phase of H i (Lall ement et al.ll200 3l). 
Recently, iHunter et al.l l)200ffl presented new ultraviolet 
Nai observations towards 74 O- and B-type stars in the 
galactic disk. Where possible they also derived column 
densities of the Nai D-lines. Then they compared the 
correlation between the Nai and Hi colur nn densities. 
While there was an excellent agreement with lFerlet et al.l 
l|l985ft for the D-lines, this was not the case for their ul- 
traviolet absorption lines. iHunter et al.l l)2006() found a 
significant offset for the correlation they got from the 
Nai D-lines compared to the relation derived from the 
UV transitions, even after accounting for saturation ef- 
fects in the D-lines. They argue that this offset is due to 
an underestimation of approximately 30 % of N(H) in the 
case of Nai D-lines observations. Recapitulating, one has 
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Fig. 5 N(H) deviations at different distances 

Shown are the differences between the model N(H) at the 
distances of the open cluster sample and the N(H[) obtained 
directly from the reddening measurements of this sample. The 
results of the extinction model based onjlakkila ct al. ( 1997) 
are plotted as grey crosses, those of the analytical model as 
black triangles. While the scatter at low distances is relatively 
small, the analytical model has a larger spread at larger dis- 
tances. 

to be aware of the still not completely solved problems 
concerning the conversion of Nai to H i column densities 
when it comes to interpreting individu al results . How- 



Fig. 6 N(H) deviations at different galactic longi- 
tudes The symbols for the models are the same as in Fig. 13 
For both models a larger scatter is visible at GLON= 0°,120° , 
and 300° 

gen as estimated bylBochkarevI (Il992l). The Local Bubble 
was taken into account by Popov et al. (2000) as a sphere 
of 140 pc radius having a constant low density of 0.1 par- 
ticles cm~'^. As noted above, there are other more recent 
theoretical 3D extinctio n models. However , the analyti- 
cal H density model bv lPopov et alJ l)2000j) seems to be 
at least as good as e.g. the underlying model used by 



ever, we regard the measurements bv iLallement et all lAmores fc Lepind l|2005j). The model of iDrimmel et all 

and the c onversion to a hydrogen column density 
according to .Ferlet et alJ l)l985j) as the best possibility at 
hand to take into account the local inhomogenity in the 
ISM for X-ray astronomy with low spectral resolution. 

At larger distances we consider two different mod- 
els - one is based completely on extinction measure- 
ments and another one is described by analytical for- 
mulae. Hakkila et al. (1992) put several extinction stud- 
ies carefully together in an easily accessible routine. All 
studies have been modified to statistically account for 
unsampled regions. Additionally, a correction method for 
the systematic underestimation between 1 and 5 kpc was 
developed. Errors were provided individually for each 
considered main survey and for their mean. The often 
large mean error values illustrate the disagreement among 

individual observations. The model by Hakkila et al. (1997J2.1 Testing the models and further improvements 
is a large scale model, capable of identifying e.g. molec- 
ular clouds at intermediate distances, but not sensitve 
to extinction variations of less than 1 degree. It is well 
suited for statistical studies 



l)2003|) cannot be applied to distances less than a few 
h undred parsecs in the S olar neighbourhood and the data 
of lMarshall et all i200(t\ has a coarse sampling with dis- 
tance bins of 1 kpc. 

Both described models are applied only for distances 
larger than 230 pc in case they provide an N(H) larger 
than derived from the sodium cube at 230 pc. Otherwise 
the hydrogen column density at 230 pc is taken since 
the column density can only increase with distance. We 
use spherical coordinates - the galactic cordinates I, b 
and distance d. The nominal sampling is one degree in 
I and b, and 10 pc for the distance. This is technically 
motivated and does not represent the actually achieved 
accuracy. Distances are covered up to 4500 pc. 



To test our models we first consider a relatively large 
number of test objects with good distances and extinc- 
tion measurements, not necessarily determined from X- 



The an alytical model we apply is descri bed bvlPopqv et alJ rays. Then we proceed to a few of the rare neutron stars 



It is based on th e formulae bv ifZ ane et alJll995t 
iDickev fc Lo ckman'l99d': lde Boeill99l|) . lPopov et alJl|200 

included additionally a galactocentric radius dependency 
for the number density of atomic and molecular hydro- 



with well known distances, having also small error bars 
;or the absorbing N{H) derived by X-ray observations. 
As iHakkila et al.l l(l997 'l note, different studies do not 
agree precisely with each other due to the various applied 
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Table 1 Testing the models for individual objects 

For each object the known distance and hydrogen column density N(H) with references are given. Furthermore are listed for 
each model: the obtained minimal distance Dmin for the lowest N(H) value, the maximal distance Dmax for the largest 
N(H) value and the expected N(H) for the known distance. The models are the same up to a distance of 230 pc (see text for 
details). The first 4 objects are neutron stars with VLBI or HST parallaxes and X-ray-obtained N(H). TY CrA lies within 
t he CrA star forming mol ecular cloud core close to R X J1856. 5-37541 HP 181 90 is thought to be situ ated behind MBM 12 . 



l: IGolden et ^llftoflh: 2: IMarshaU fc Sdigi3 112001; 3: IDodson et "!!! pno^) : 4:IPavlov et alJ (l200l : 5 :ICaraveo et alJ 
6: IJackson fc Halpcrn 1120051) : 7: IKaplan et al.l (120021k 8: our EPIC-pn fit; 9: ICasev et al.l (119981k 10: IForbrich etalT 
ll: IHobbs et al. (,.198&) : 

" Zog(N(Nal)) = 10.6 < 11 IFerlet et~an ([T985ll formula was doubted bv IWeltv eFaLl EoM) . see text; '' note the recently 
refined parallax bv lvan Kerkwiik fc Kaplan! II2006I) yielding a distance of I6III4 pc; mean value; lower limit for N(Nal) 
was converted to N(H); out: needed N(H) is not reached up to the here considered 1 kpc distance 



Name 


Distance 


Ref. 


N(H) 


Ref. 




analytical model 




extinction model 














N(H) 


Dmin Da/ax 


N(H) 




[pc] 








[pc] 


[pc] 




[pc] 


[pc] 




PSR B0656+14 


288lp 


1 


1.73 ±0.18 


2 


240 


260 


2.68 


230 


240 


2.86 


Vela 


287+^^ 


3 


3.3 ±0.3 


4 


300 


330 


2.73 


860 


910 


1.07 


Geminga 


15711? 


5 


1.76 ±0.95 


6 


220 


300 


0.09" 


220 


240 


0.09" 


RX J1856.5-37541 


140 ± 40'' 


7 


0.74 ±0.10 


8 


130 


140 


0.87 


130 


140 


0.87 


TY (JrA 


129 ± 11 


9 


130 ± 0.3" 


10 


out 


out 


0.08 


out 


out 


0.08 


HD 18190 


185 


11 


> 7.21^* 


11 


440 




5.05 


800 




5.05 



methods or objects studied. The latter are usually very 
coarsely scattered and measurement errors can influence 
the extinction values of a large region. One possibility 
to overcome this problem is investigating open clusters. 
An open cluster has the advantage that one can measure 
a number of stars at approximately the same distance, 
which is usually well known. Therefore, open clusters 
are a good choice in aiming to minimize the extinction 
measurement errors along one line of sight. Usually one 
measures the reddening E{B — V). We concentrate here 
on the recent compr ehensive reddening mea surements for 
650 open clusters bv lPiskunov et al.l l|200^. representing 
a complete sample up to distances of about 850 pc from 
the Sun Details about this open cluster study can be 



9 • 10 cm~ , the corresponding mean value for the ana- 
lytical model is 27-10^° cm~^. Considering only distances 
up to 1 kpc these values are 3.3 ■ 10^° cm~^ and 8.1 • 10^° 
cm~^ respectively. Interestingly at some galactic longi- 
tudes {I « 0°, 120°, and 300°; see Fig. and at the 
galactic latitude 6 « 0° t he scatter is more pro nounced. 
The extinction routine bv'H akkila et alJ ()l997j) predicts 
the observed properties of the open clusters better than 
the analytical model. This can be partl y explained by the 
use of (older) op en cluster data (e.g. iFitzgeraldl Il968() 
in the routine by iHakkila et al.l ()l997(l . The analvtical 
model tends to overpredict the extinction by an order of 
magnitude, especially at larger distances. Due to the con- 
vincing advantages of extinction data derived from open 



found in lPiskunov et all l)200(J ); lKharchenko et alJ ()2005alblfc lusters. we include them finally as local inhomogenities 
To conver t the reddening E(B — V) we apply the formula in our models for distances larger than 230 pc. Again, 
bv lPares ce (198^: N{H) = 5.5 x 1Q^^E{B - V) [cm'^j . we take into account that the column density can only 
This is quite similar to the slope found in the X-ray study increase with distance. If N(H) outside 230 pc is lower 
for the extinction At/ bvlPrcdchl fc Schmitt ( 1995): N(H) ^han the value at 230 pc in this dire ction we conse- 

uently changed the values based on the lLallement et all 
to the lower value until an N(H) increasing homo- 



1.79±0.03xl02iAy-0.41 [cm-2J if Ay = 3.1E{B-V). 
This is not ed here because there is a d eviation of the cor- 
relation bv lPredehl &: Schmittl l)l995() between Ay, indi- 
cator of the dust, and N{H), indicator for cold gas and 
dust, at low distances. This is probably due to a signifi- 
cantly lo wer amount of dust in the Local Bubble, as also 
found bv lPredehl fc Schmitti (|l995l) for the source LMC 
X-1. Therefore, we consider only open clusters with dis- 
tances larger than 230 pc. 

We compared the hydrogen column densities we inferred 
from the open cluster reddenings with those o btained by 



geneously with distance is reached. 
Aiming to apply our N(II)-models to X-ray detected neu- 
tron stars their performance for well known sources is 
interesting. We selected four of the rare NSs having at 
the same time known parallaxes and hydrogen column 
densities with small error bars obtained from X-ray ob- 
servations. In Tabled we present the model-derived dis- 
tances for the measured N(H) range as well as the model- 
derived N(H) at the given parallactic distance. While 



the mo dels. 628 out of the 650 open clusters bv lPiskunov et tiik analytical model gives reasonable results for PSR 
(|2006l) he within the considered data cube with distances B0656+14, Vela and RX J1856. 5-37541 ( RXJ1856 in the 
larger than 230 pc and smaller than 4.5 kpc. The scatter following), the extinction model based on lHakkila et alj 
in the obtained extinction differences increases with dis- l|l997l) underestimates the column density towards Vela 
ta nces for both models ( see Fig.jS}. For the model based significantly. As both ISM models are equal up to 230 pc, 
on IHakkila, et all lll997li the mean deviation is around both give the same minimal distance of 220 pc towards 
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Geminga whic h deviates from the claimed distance by 
ICaraveo etahl (|I996). Interestingly Walter et al. 2006 
(this volume) reported recently a revised distance of 
254^gg^ pc towards Geminga. In the case of RX J 1856 the 
derived minimal and maximal distances are completely 
includ ed in the N(II) cube obtained by the lLallement et all 
l)2003(l measurements. Here, an error of 25 pc applies 
while an error estimate is difficult for the analytical model 
at distances larger than 230 pc. The Hakkil a et al. (J^997} 
model has often a formal mean error of the same order 
as the obtained N(H) values. 

As noted in above it is still debated how good Nai 
D-lines are as tracer for hydrogen. In particular it is un- 
clear whether molecular hydrogen is well traced. In Ta- 
ble ^ we include the results obtained for two further test 
objects situated in or close to molecular cloud cores - 
the star-forming site C rA close to RXJ1856 and MBM12 
(|Magnaniet alJll985^ . While towards MBM12 there is 
an enlarged but still underestimated amount of hydrogen 
(see also ,Lallement et al. 2 003 for discussion of MBM12), 
the small CrA cloud is clearly missed. Both examples in- 
dicate that one has to be cautious i n applying the N{H) 
cube obtained from the lLallement et al. (20(^ measure- 
ments in the direction of small molecular clouds. How- 
ever, there are only few such density enhancements of 
molecular hydrogen in the close Solar neighbourhood. 
The extinction model as well as the analytical model are 
in principle sensitive to molecular hydrogen. At large dis- 
tances > 1 kpc the extinction model may miss molecular 
clouds due to their low angular size and bad sampling 
while the analytical model does not consider individual 
clouds. 



Table 2 Distances obtained for the M7 

The N(H) was obtained by blackbody and absorption line 
(tabs) XSPEC-&ts from XMM-Newton EPIC-pn observa- 
tions (abundances and cross-sections are chosen according to 
^ilms ct al. 2000). The number of lines used is indicated in 
parenthesis. The corresponding model prediction of the dis- 
tance is given for the analytical model: dana, the extinction 
model based on Hak kila et al.l (|l99t) dext, and the analytical 
model plus consideration of possible 30% underprediction by 
the sodium D-lines: danaiso- See also text for discussion. 



3 Applications of the models 

3.1 Distance estimates 

A possible application is to estimate the distance to neu- 
tron stars without parallaxes or other distance measure- 
ments. This is the case of the majority of the M7 where 
despite for RXJ1856 only a preliminary parallactic dis- 
tance of 330^1™ pc is reporte d for RX J0720.4-3125 by 
Ivan Kerkwiik fc KaDlanll2006l (this volume). We obtain 
N(H) by considering the best XSPEC-fits for black- 
body and absorption lines of all currently available XMM 
EPIC-pn observations reduced with XMM SAS 6.5 for 
each neutron star. As noted by Haberl 2006 (this vol- 
ume), in this respect the XMM EPIC-pn is the best- 
suited instrument for relative measurements with rea- 
sonably small errors (~ 0.1 • 10^" cm~^ for all objects). 
The model predictions for these column densities can be 
found in Table [5] As noted before, the technical sam- 
pling is 10 pc. If a N(H) value lies between two sam- 
pling points we indicate this by noting e.g. 235. How- 
ever, even errors in the best-studied regions (< 230 pc) 
are around 25 pc , and can be much larger otherwise. 
In cas e of the extinction model based on Hakkila et al. 
l)l997|) one cannot always derive a convincing value due 



name 


N(H) (#lines) 




dext 


danalSO 




[lO'o cm-2] 


[pc] 


[pc] 


[pc] 


RX J1856.5-3754 


0.7 (OL) 


135 


135 


125 


RX J0420.0-5022 


1.6 (IL) 


345 




325 


RX J0720.4-3125 


1.2 (IL) 


270 


235 


265 


RX J0806.4-4123 


1.0 (IL) 


250 


235 


240 


RBS 1223 


4.3 (IL) 








RX J1605.3+3249 


2.0 (3L) 


390 




325 


RBS 1774 


2.4 (IL) 


430 




390 



to bad sampling with the (low) extinction not changing 
over a large scale of distances (e.g. up to 1 kpc). For 
the analytical model only the high latitude object RBS 
1223 is a problem where at 1 kpc one arrives only at 



3.2 • 10 cm~ . The limitations of this large-scale model 
appear here. However, absorption lines in the spectrum 
also influence N(H). Schwope 2006 (this volume) fitted 
the spectrum of RBS 1223 using recent new observations 
and obtained N(H)= 3.7 • 10^" cm~^ for one absorption 
line. This value is comparable to our result in Table El 
taking into account the additional data. However, the fit 
was not very good and Schwope 2006 obtained better fits 
with two absorption lines N(H), yielding values ranging 
from 1.2 • lO^*^ cm^^ ^ g . ;lo20 cm-^. The lower value 
would correspond to a distance of 525 pc when using the 
analytical model. 

Overall, the distances derived from both models towards 
two neutron stars are in acceptable agreement: 
for RX J0720.4-3125 with 250 pc and for RX J0806.4- 
4123 with 240 pc. Since both NSs are close to the bor- 
der of the sodium measurements, an error of « 25 pc 
seems reasonable. In general, the values of the analytical 
model have to be used with caution as local dumpiness 
of the ISM is not considered by this model at distances 
> 230 pc. Considering the open cluster test from above, 
the mean N(H) deviation of the analytical model up to 
500 pc is 0.05 • 10^° cm~^; the standard deviation, how- 
ever, is 4.9 • 10^*^ cm~^. 

3.2 Population synthesis 

The main motivation for us to take into account bet- 
ter absorption treatment is an improved population syn- 
thes is model based on th at of Podov et al. ( 2000, 2003) 
and iPopov et alJ l|2005c|) . A detailed discussion of the 
improvements is beyond the scope of this article and 
will be provided elsewhere (Popov, Posselt, in prepara- 
tion). Therefore, we summarize here oidy the new devel- 
opments and concentrate on illustrating the importance 
of the interstellar absorption. 
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Fig. 7 The logA*"— logS curves for different models 

Each curve shows the expected number of observable neutron 
stars at ROSAT PSPC countrates. The solid line represents 
the curve for the new progenitor distribution but old absorp- 
tion model; the dotted and the dashed line are the results 
of the old simulations by e.g. f'ouov et al. (2003) for Gould 
Belt radii of Rgb ~ 500 pc and Rob = 300 pc respectively. 
The only difference to the solid curve is the progenitor distri- 
bution. The results of taking into account the new progeni- 
tor distribution and additionnally the new absorption models 
are shown by the dashed-dotted line for the refined analytical 
ISM model and the dot-dot-do t -dash ed line for the extinction 
model based on lHakkila et sL (1997). Furthermore the mea- 
surement points of young N Ss with thermal X-ray emission 
are plotted as in lPopov et al.l Ii2003ti . For discussion see text. 

While iPopov et all l)2003l l2005cj) modeled the neutron 
star progenitor distribution as coming from infinitely 
thin disks, either from the galactic plane or the Gould 
Belt, the new initial distributio n is more r ealistic. Up 
to the Hipparcos limit of 400 pc l|ESAlll997() the known 
B2-08 stars are considered individuall y as well as their 
affiliation to an OB association (d c Zeeuw et al. 1999) 
resulting in birth properties depending on the age of 
the OB-association. The Gould Belt birth rate is applied 
here. For distances above 400 pc, the galactic disk is con- 
sidered f or few neutron stars randomly as well as 36 asso- 
ciations l|Blaha fc Humphrevslll989t iMel'nik fc Efremovl 
for most of the neutron stars. Due to the un- 
known ages of the associations, the birth probability 
is set proportional to the number of association mem- 
bers. In Fig. [TJ the logAi'-logS' curves are shown for the 
old thin-disk Gould Belt models with cither 300 pc or 
500 pc radius, the new progenitor distribution as well 
as for the the new progenitor distribution considering 
both absorption models. All si mulations were don e for 
the same mass spectrum as in iPopov et al.l ()2006j) , the 
cooling curve set labeled Model I in the same paper, the 
sa me abundances and cross -sections (in this case those 
of iMorrison fc McCammonI (^983) for comparison) and 
identical remaining parameters. The results for the new 
progenitor distribution lie in-between those of the two 
Gould Belt sizes, a little apart from the measured inte- 



grated number of objects. This is also true for the curve 
considering the absorption by the analytical model. At 
the bright end the differences introduced by the ana- 
lytical absorption model are only very small. Towards 
lower countrates, the curve lies significantly below the 
one without the refined 3D ISM model. Considerin g the 
absorp tion based on the extinction model by Hakkil a et alJ 
l)l997f) gives another picture. At the bright end the curve 
is much lower than that without the refined 3D ISM 
model, while the difference becomes smaller towards low 
countrates. We note here again that both absorption 
models a re the same up to 230 p c, based on the measure- 
ments bv lLallement et all l|2003() . Both also consider the 
open cluster data, however slightly differently as the col- 
umn density is only allowed to increase. The curve of the 
analytical model shows best the influence of an improved 
ISM distributio n model cornpared to the old absorption 
treatment (e.g. IPodov et al. | l2003() where th e analy tical 
model without the data bv'Lall ement et al.l l)2003() and 
Piskunov et al. ( 2006) was used) . One might have ex- 
pected that differences should be largest at low distances 
and hence for bright sources. However, averaging over 
the small Local Bubble gives apparently the same result 
as the simpler description before. The influence of the 
better model for the Solar vicinity is apparent in the ex- 
pected location likelihoods for neutron stars. Going to 
larger distances the additive column densities are espe- 
cially large around I « 300° to / w 60° in the galactic 
plane where many objects could be expected. Thus, the 
number of predicted neutron stars is reduced. The N(II) 
model based on the extinction study by Hakkila et al.l 
II1997D shows on average much higher N(H) already at 
240 pc compared to the analytical model. Given the 
nearly rectangular shape of the regions with relatively 
high values in the I, b plane these seem to be caused at 
least partly by bad sampling and the necessary statis- 
tical treatment. The higher N(H) results in less observ- 
able sources bringing the corresponding logA^-logS* much 
closer to the actually measured number of sources at the 
bright end. 

We made similar test plots for the expected neutron star 
nu mber per square degre e like the one presented in Fig. 6 
bv iPopov etal] l|2005c!). We summarize here only the 
main results, for detailed discussion see Popov, Posselt 
2006 {in prep.). Contrary to the old model the predicted 
neutron star number is enhanced in regions including the 
nearby Gould Belt OB associations (e.g. Upp Sco) with 
the new progenitor distribution. Sources are not any- 
more homogeneously distributed at all longitudes along 
the planes of the galactic disk or the Gould Belt. When 
it comes to consider the 3D ISM models, the influence 
of the absorption is on much smaller scales. Less sources 
would be expected at higher latitudes, and towards some 
small regions the high extinction hinders the possible 
detection of neutron stars. Differences between the two 
absorption models are pronounced only in small regions 
which to discuss in detail is beyond our intention here. 
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4 Conclusions 



An understanding of the interstellar absorption is cru- 
cial for the interpretation of low-resolution soft X-ray 
spectra. Taking into account the 3D distribution of the 
ISM is important to determine the absorption of known 
sources as well as for searches of new soft X-ray sources 
towards a particular direction in the sky. We presented 
two simple models for 3D ISM distributions together 
with possible applications: distance estimation and the 
population synthesis for young isolated thermal, thus soft 
X-rays emitting neutron stars. Both models are large 
scale models and have to be used with caution especially 
considering distance estimations. Concerning population 
synthesis the two models influence the logNAogS curve 
discriminatively, while the differences for the predicted 
spatial distribution of observable neutron stars are only 
small. 
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